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Abstract 

We show that it is now possible to image optically thick Lya clouds in fluores- 
cent Lya emission with a relatively long (~ 20 hr) integration on a large (~ fOm) 
telescope. For a broad range of column densities {N ^ 10"^^'^ cm~^), the flux of 
Lya photons from recombination cascades is equal to ~ 0.6 times the flux of ion- 
izing photons, independent of the geometry of the cloud. Additional Lya photons 
are produced by coUisional excitations when these are the cloud's primary cool- 
ing mechanism. For typical physical conditions expected in optically thick clouds, 
these mechanisms together lead to a Lya emission flux that is ~ (2/3) {p) j vq 
times the flux of ionizing photons, where (^v) is the mean frequency of ionizing 
background photons and vq is the Lyman limit frequency. Hence, measurement 
of the surface brightness from an optically thick cloud (known to exist, e.g., from 
a quasar absorption line) gives a direct measure of the energy in the ionizing ra- 
diation background. Moreover, in the same long slit spectrum one could hope 
to detect emission from ~ 200 other Lya systems. Such detections would allow 
one to make a 2-dimensional map of the Lya forest. By taking a series of such 
spectra, one could map the forest in three dimensions, revealing structure in the 
high-redshift universe. 
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1. Introduction 



Resonant absorption against the continuum of a bright background source is an 
exquisitely sensitive way to detect diffuse gas. The Lya forest in the spectra of high 
redshift quasars provides a detailed picture of structure in the universe at 2; ~ 2 — 5, 
revealing systems that are much more common than galaxies but whose internal 
gas densities are far lower. Each quasar spectrum yields a f-dimensional map of 
the neutral hydrogen distribution along its line of sight. Our limited information 
about the sizes and morphologies of Lya absorbers comes from the rare instances 
where two nearby quasars (or the multiple images of a gravitational lens) probe 
neighboring lines of sight. In this paper, we argue that lO-m class telescopes make it 
feasible, though still challenging, to detect the optically thick end of the Lya forest 
(neutral column density N ^ IO"^^cm~^) in Lya emission, thereby building up a 
3-dimensional picture of the properties and space distribution of these systems. 

There have been many attempts to detect Lya emission from high redshift 
"primeval galaxies." Some of these efforts have covered large areas of sky (e.g., 
Thompson, Djorgovski, & Trauger 1995; Thompson & Djorgovski 1995), while 
others have targeted known damped Lya absorbers (e.g., Lowenthal et al. 1995), 
which are plausible progenitors of disk galaxies. The motivating idea behind these 
searches is that regions of active star formation should produce large numbers of 
Lya photons. The paucity of successful detections has therefore been interpreted 
as evidence for either low rates of star formation or internal dust extinction in these 
systems. However, hydrogen clouds should emit some Lya radiation even if they 
have no ongoing star formation at all, because the metagalactic UV background 
photoionizes the clouds (or at least their optically thin skins), and many of the 
resulting recombinations produce Lya photons that ultimately escape. 

The idea of searching for this fluorescent emission from Lya systems has been 
discussed previously, e.g. by Hogan & Weymann (1987) and by Lowenthal et al. 
(1990), who used their null result to place an upper limit on the UV background 
intensity. Hogan & Weymann (1987) estimated the surface brightness of optically 
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thin {N ^ 10"^ cm~^) clouds and concluded that these might plausibly be detected 
with 2-m class telescopes. Their conclusions were substantially more optimistic 
than those reported below, largely because they assumed a fiducial UV background 
intensity more than 100 times higher than the more recent observational estimates 
that we adopt. With the lower background intensity, optically thin clouds are 
difficult to detect even with 10-m telescopes. We therefore emphasize the possibil- 
ity of detecting optically thick systems, which, roughly speaking, absorb the UV 
continuum radiation and mirror it back into space in the concentrated form of a 
double-peaked Lya line. We discuss a number of physical effects that infiuence the 
optimal design of search strategies and the interpretation of detections and null 
results. 

Ionizing radiation gives rise to Lya emission in two distinct ways. First, in 
ionization equilibrium, each ionization leads to a recombination, and most of these 
in turn ultimately yield a Lya photon. Second, in thermal equilibrium, the kinetic 
energy of the ionized electron heats the gas, and this heat is ultimately emitted in 
radiation. Most of this radiation is in the form of free-free emission, free-bound 
emission, or line emission induced by coUisional excitation. At high densities and 
temperatures r~2 — 5xlO^K, conditions that are likely to prevail in optically 
thick clouds, line emission is the dominant cooling mechanism, so a large fraction 
of the deposited thermal energy also comes out in Lya photons. 

Our approach in this paper is first to analyze the recombination Lya photon 
production, which depends primarily on the total number of ionizing photons. 
We give a simple analytic treatment in § 2, then present the results of numerical 
simulations in § 3, where we focus on the frequency and angular distribution of 
the emitted line photons. In § 4 we analyze the effects of cloud topology (coherent 
absorbers vs. clusters of cloudlets) on the expected form of the emission signal. 
In § 5, we turn to the excitation Lya photon production, which depends primarily 
on the total energy of the ionizing photons. We show that, for optically thick 
clouds, the total Lya emission is proportional to the total energy of the ionizing 
background. Finally, in § 6, we present a signal-to-noise estimate for a fiducial 
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observation and discuss what could be learned from various types of searches for 
emission from Lya systems. 



2. Recombination Emission from Lyman a Clouds 



Consider first a Lya cloud which is optically thin to ionizing radiation, N <^ 
10"^^cm~^. The hydrogen atoms will then be ionized at a rate F, 



where J(z^) is the photon fiux density (in ergcm"^ s~"^ sr~"^ Hz~"^), cr(z/) is the ion- 
ization cross section, and hpo = 13.6 eV is the Rydberg energy. We assume that 
the cloud is in ionization equilibrium, i.e., the ionization rate equals the recom- 
bination rate, which should be a good approximation for Lya clouds. A typical 
recombination leads ultimately to the production of a single Lya photon, and this 
photon resonantly scatters off H atoms until it escapes the cloud. There are, how- 
ever, several loss mechanisms. First, for a broad range of relevant temperatures 
centered on T = 10,000 K, ~ 38% of recombinations go directly to the ground 
state (Osterbrock 1989). These recombination photons escape an optically thin 
cloud without generating Fya photons. Second, a fraction ~ 32% (which is again 
only weakly dependent on T) of recombinations ultimately produce excited H in 
the 25* rather than 2P state (Osterbrock 1989). These atoms decay by 2-photon 
rather than Fya emission. Hence a fraction ?7thin ~ 42% of recombinations ul- 
timately yield a Fya photon. If the cloud h a cosmologically nearby 
observer [z <^ 1) at distance D would see a fiux of Fya photons ?ythin^^^r/47r_D^, 
corresponding to a number fiux (photons cm~^ s~"^ sr~"^) of $obs = '?thin^r/47r. At 
cosmological distances, bolometric surface brightness is reduced by (1 -|- z)~'^, so 
photon number density is reduced by {1 -\- z)~^ . Hence the observable Fya photon 
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flux is, 

$obs = (1 + ^)"SthiniV / diy^{iy)a{iy) (optically thin), (2.2) 



where ?7thin ~ 42%. The cross section for hydrogen ionization is (Osterbrock 1989), 



1 — exp(— 27r/e) 
where 



z/Q / tan e 

— exp 1 

1^ \ e 



4 /X -2.75 

, (2.3) 



/ \ 1/2 

(70 = 6.3 X IQ-^^cm^ e = f — - 1 j . (2.4) 

Consider next a spherical cloud of radius r which is optically thick to ionizing 
radiation. Every incident photon with p > will ionize a hydrogen atom, which 
will then recombine and, in most cases, produce a Lya photon. It is still true 
that ~ 38% are captured directly into the ground state, but most of the resulting 
continuum photons are captured again. We flnd numerically in § 3 that only ~ 5% 
leave the cloud without converting into a line photon. As in the optically thin case, 
~ 32% of line photons end in 2-photon decays. In § 3, we estimate the fraction of 
the remaining Lya photons which are absorbed by dust before they escape. For our 
flducial cloud with N ~ 10"^^ cm~^, this fraction is ~ 4%. The cloud will therefore 
emit Lya photons at a rate ?7thick ~ 62% times the rate at which it absorbs ionizing 
radiation, i.e., 47r^?ythick'"^ / dvcl>[v). (The parameter ?7thick has a weak dependence 
on HI column density, as shown in § 3, below.) We now assume that the cloud 
has a "sharp" boundary, so that the region of continuum absorption is the same as 
the region of Lya emission. Making use of this rather strong assumption, whose 
validity we discuss in § 3, we repeat the above calculation and flnd. 



^obs = (l + ^;) Nthick j dv(j){v) (optically thick). (2.5) 

While this calculation was performed for a spherical cloud, it actually remains valid 
for any sharp-edged optically thick cloud. This conclusion follows from detailed 
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balance and the fact that both the ambient radiation and the emission process are 
(by assumption) isotropic. Note that we have ignored ionization of Hell by ener- 
getic photons [p > 4z/o)- We include Hell ionizations in our numerical simulations 
reported below, but their net effect is small. 

3. Lya Scattering 

In an optically thick cloud, ionizing photons will be absorbed at an (ionizing) 
optical depth ~ f , corresponding to column depth N ~ [p /po)^'^^crQ^ ~ 10"^^ cm~^, 
where we have taken p / ~ 2. This corresponds to an optical depth r ~ 10^ for 
the emitted Lya photons, if the cloud's I-dimensional, rest-frame velocity disper- 
sion (which may reflect both thermal broadening and bulk internal motions) is 
a ~ 35kms~"^. Scattering of Lya photons in an optically thick hydrogen cloud 
has been studied for at least half a century (Holstein 1947; Neufeld I99I and refer- 
ences therein). However, the geometry of the present problem differs significantly 
from that considered in most previous treatments. In particular, our Lya photons 
are generated in the cloud's "shielding" layer, and since this is likely to be small 
compared to the cloud, the appropriate geometry is a slab. We therefore focus on 
the problem of the distribution in redshift and angle of photons which escape from 
such a slab. The bottom line is rather simple: a typical Lya photons "scatters" 
(more precisely, it is absorbed and reemitted) until it is emitted by an atom with 
velocity ~ it y^2 ln{T)a ~ ±4(7 relative to the central velocity of the cloud. The 
emission is isotropic, and the profile is double peaked at the ~ 4 cr wings of the 
cloud. 

In the physical conditions of Lya clouds, the time between photon absorption 
and emission is very short compared to the collision time and can be neglected. 
Hence, the absorption and emission together can be regarded as scattering. The 
scattering is not quite isotropic because of the sin^ ^ dependence of the 1S-2P 
transition, but this anisotropy has virtually no effect. As we show below, the 
typical column densities at escape are N ~ 10"^^ cm^, so the damping wings of the 
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Lya line are unimportant. Hence, the photon is always absorbed by an atom with 
the same velocity projected along its direction of motion as that of the atom which 
emitted it. 

Although the escape of Lya photons from a cloud is often called "diffusion", 
this process is not primarily diffusive in either physical or frequency space. Rather, 
the photons tend to maintain their initial spatial distribution, and up to one or two 
scatters before a photon finally escapes, its redshift is near the mean of the cloud. 
Each photon in this distribution suffers a single catastrophic scatter which suddenly 
moves its dimensionless Doppler shift A to the extreme tail of the distribution, 
A ~ ±4. Here, 



where 6p is the Doppler shift, p is the line center, and a is the f -dimensional 
dispersion of the gas. The photon then faces a very low effective optical depth 
Tgff = Texp(— A^/2) where r is the optical depth at line center. It therefore 
escapes directly from high r. This simple model of the escape process, justified 
analytically and numerically below, allows one to make an accurate estimate of the 
distribution of emitted photons as a function of redshift and angle. In brief, for 
an initial distribution at characteristic optical depth tq, the photons are emitted 
with angular dependence oc cos 6 dcos 6 (just as is the case for diffusive escape) and 
with redshift sharply peaked at it Apeak? where Apeak ~ y^2 ln(To sec 9) (which is 
substantially different from diffusive escape). The width of each of these peaks is 
~ Apg^]^. Here 9 is the angle of escape relative to the normal of the surface of the 
cloud, which is assumed to have little curvature on scales of the physical depth at 
which the photons are captured. 

After any given scatter, the probability that the photon is shifted by A is 
p{A) = (27r)~"^/^e(A), where e(A) = exp(— A^/2). Such photons will encounter a 
reduced density of hydrogen scatterers and will therefore travel a factor farther 
than photons at the line center before colliding again. [The density of wing photons 
in the gas at any given time is therefore uniform in A, i.e., with no Boltzmann 
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factor, exp(— A^/2).] If one of these photons is headed in the direction ^, then it 
will face an effective optical depth Tgff = er sec 6 and so will escape with probability 
exp(— Teff). Hence, at fixed optical depth r, the rate of escape F is, 

fir 

oc ee-^^^^^^ (3.2) 



dcos 9dA 

At fixed ^, the peak of this distribution is at 



Apeak(^) = V21n(Tsec^), (3.3) 
and the integrated emission rate, 

dV cos 6 

dcos 9 T Apeak ' 

is very nearly proportional to cos 6. The width of each of the two peaks is ~ A~J^]^. 

The distribution of optical depths of the Lya photons prior to escape can 
be roughly characterized as an exponential, exp(— t/to)c/(t/to). We focus on the 
emission at fixed angle, 6. Since the distribution in A at fixed r is comprised of 
two narrow peaks, we approximate it as two (5-functions at A = ibApeak- From 
equation (3.3) and the definition e = exp(— A^/2) we find that r = cos 9, and 
this substitution leads to a total emission rate. 



"^^'"^ oc Ae-ie-^°^^/^^". (3.5) 



dcos 9dA 



This distribution is narrowly peaked around A ~ ibY'21n(To sec^). The net result 
is that ionizing photons penetrate to a typical (Fya) optical depth r ~ 10^, then 
the emission takes place from a similar depth, with an angular distribution very 
similar to that from a diffusive optically thick surface, and with a shift distribution 
composed of two sharp peaks at A ~ ±4. 
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We confirm this picture by means of a Monte Cario simuiation. We first ir- 
radiate an optically thick slab of hydrogen {N = iO"^^cm~^) with an ionizing 
spectrum taken from Haardt & Madau (i995; hereafter HM). HM compute the 
ambient UV background that would be produced by the emission from observed 
quasars after it is filtered through the Lya forest; for the calculations in this pa- 
per, we use HM's z = 2 spectrum. (We assume that for p > 4z/o, the ionization 
cross section is 12.5(Th(?^/4), corresponding to a ratio of Hell to HI of 50. How- 
ever, this assumption has relatively little effect because the HM spectrum has very 
few such photons.) We determine the fraction of photons which pass through the 
cloud without interacting (0.5%) and the fraction which generate Lyman contin- 
uum photons that escape without forming line photons (5%). Next, we allow each 
resulting Lya photon to scatter until it leaves the cloud, and we tabulate A^walk? 
the total column density through which the photon randomly walks in the course of 
escaping. We then compute the probability that each photon is absorbed by dust, 
1 — exp(— A^^alk/^*)? where = 10^^ cm~^. This value is adopted based on the 
work of Fall & Pei (1995), who argue that damped Lya systems of column density 

typically extinguish about 1 mag of background quasar light. We find that 
dust absorption is small, ~ 4%. Of course, the absorption could be still lower if 
the dust-to-gas ratio for N ~ 10"^^ cm~^ systems is lower than for damped systems. 

Figure 1 shows the cumulative angular distribution of the emitted photons com- 
pared to the curve (cos^ 9) one would expect for optically thick diffusive emission. 
There is virtually no difference. 

Figure 2 shows the dimensionless Doppler shift A summed over all emission 
angles. As expected there are two narrow peaks. We have confirmed (but not shown 
in the Figure) that the distributions at fixed emission angle have even narrower 
peaks, whose separations grow oc Y^ln(T sec 9). 

Figure 3 illustrates the essentially non-diffusive nature of the emission process. 
The solid and solid bold histograms show the initial and pre-emission distributions 
of the photons in optical depth. The two distributions are broadly similar. It is true 
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Figure 1. Cumulative distribution of emitted Lya photons as a function of the cos of the angle 6 relative 
to the normal of the surface (solid curve). The emitting cloud is assumed to have a slab geometry with 
a column N = lO^^cm"^. Shown for comparison is the distribution expected for optically thick diffusive 
emission (dashed curve). The fact that the two curves are almost identical means that the observed intensity 
of Lya photons does not depend on the orientation of the cloud surface. 



that the photons are typically emitted at about half the optical depth at which they 
are absorbed, but since the characteristics of emission are determined by Y^ln(T), 
this reduction does not qualitatively alter our basic picture. The small role played 
by diffusion can be understood as follows. The distribution of Doppler shifts is 
~ e(A). Typically a photon will directly escape after N ^ t collisions by being 
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A = (Av/v)/(a/c) 



Figure 2. Distribution of photons as a function dimensionless Doppler shift, A, of the Lya photons emitted 
by an optically thick (N = IQ-'^^cm"^) cloud. 

scattered to Apeak ~ a/2 In r. We can estimate the number of mean free paths 
(mfps) traveled by the photons at all shifts |A| < Apeak- After each collision, 
the mean square number of mfps is 2e~^ . Consequently, the total mean square 
movement is ~ ~ T . That is, the photon diffuses ~ r mfps, 

comparable to the number it moves in its largest single scatter. The order of this 
effect is clear in Figure 3. If the escape were primarily diffusive, then most of the 
photons would escape from r ^ 1. By contrast, we find that ~ 1% of the photons 
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Figure 3. Initial (solid histogram) and pre-emission (solid bold histogram) distributions of photons in units 
of one Lya optical depth. The distributions are broadly similar. The main difference is that photons tend 
to diffuse a factor ~ 2 closer to the edge before they escape. At N = IQ-'^^cm"^, a fraction ~ 87% leave 
from the same side (left) that they entered. The remaining 13% escape from the far side, including ~ 4% 
that escape from the right-most bin which contains ~ 2 x 10'* optical depths. 



escape from regions of r < 1, ~ 5% from r < 10 and ~ 14% from r < 100. These 
values correspond to column densities of ~ 10 cm" -2, N ~ lO^^cm-2, and 
N ~ 10"^^ cm~^. Thus, even if an optically thick Lya cloud has an extended halo 
of relatively low density gas, the emission of Lya comes primarily from the high 
T region where the ionizing photons are absorbed. This validates the assumption. 
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made in § 2, that the clouds can be considered "sharp-edged". 

Unhke Lya photons created in star forming regions, the Lya photons created 
by the ionizing background are virtually unaffected by dust. This is because the 
ionizing background penetrates only to a depth N ~ f 0"^^ cm~^, so the Lya photons 
have only to traverse this depth in order to escape. By contrast, the photons 
generated by hot stars deep within a damped Lya system {N ~ fO^"^ cm~^) must 
travel through a column density ~ fOOO times larger. 

Figure 4 shows the relative importance of various loss mechanisms — passage 
through cloud, direct capture to ground state, cascade through two-photon con- 
tinuum, and extinction by dust — for slabs with column densities N = 10^'^ — 
10^^ cm~^. The highest line shows the fraction of ionizing photons which are ini- 
tially captured by the cloud. The next three lines are the fractions which produce 
line photons, which produce Lya photons, and which produce Lya photons that 
escape before hitting a dust grain. At = fO^*^ cm~^ virtually all (> 98%) of the 
photons which escape come out the side that they entered. This means that for 
all N > 10^^ cm~^, the results will be the same as for the final point. Thus, the 
emitted fiux peaks at ~ fO"^^ cm~^, but it is close to the peak for all column 
densities A^ k, fO^^'^ cm'^ 

4. Effects of Cloud Topology 

In § 3, we modeled the cloud as having a slab-like shielding area with a single 
characteristic velocity dispersion. An alternative model for optically thick clouds 
is that they are composed of many cloudlets orbiting in a halo of hot, thermally 
unstable gas (e.g. Mo 1995). The observed velocity dispersion would then be due 
to the bulk motion of the individual cloudlets rather than the internal motions. 
Thus, the observed absorption line would be a composite of absorption from several 
individual cloudlets (the halo might contain many cloudlets, but only a few would 
intersect a given line of sight). 
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Figure 4. Fraction of ionizing photons which ultimately lead to Lya emission as a function of cloud column 
density N. A slab geometry is assumed. The four curves shown represent the fraction of ionizing photons 
which lead to hydrogen ionizations, line photons, Lya photons, and emitted Lya photons. 

Introducing this model has virtually no effect on the surface brightness of Lya 
emission. It is still the case that essentially all ionizing radiation incident on the 
cloud will be absorbed. The Lya yield from this absorption will be essentially 
unchanged. What will change is the profile of emitted radiation. Consider an 
observation of emission toward a quasar with a known system along the line of 
sight. To be specific, let the individual cloudlets have dispersions a ~ lOkms""^ 
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and be moving with respect to one another at a' ~ 70kms~"^. Suppose that the 
spatial resolution is good enough to resolve the individual cloudlets in the cloud 
system. What would one see directly on the line toward the quasar? First, one 
would see emission with a profile of the form shown in Figure 2, two sharp peaks 
separated by ~ 80" ~ 80kms~"^. This would come from the foremost cloudlet 
along this line of sight. There would be other cloudlets behind it (accounting 
for the observed dispersion in the absorption line), and these would emit similar 
double-peaked profiles (with different line centers). However, if any of these peaks 
fell between the two observed peaks from the front most cloudlet, these peaks would 
be absorbed. If any of the cloudlets had peaks outside this range, they would be 
transmitted through the front most cloud. Hence the overall spectrum would have 
a gap of ~ SOkms""^ surround by two (or possibly more) peaks. 

If one could not resolve the individual cloudlet surfaces, either because of inad- 
equate angular resolution or because of poor signal to noise, then one would see a 
superposition of such structures each with a different center. If the cloudlets were 
at high Mach number (large bulk velocities compared to internal dispersion), the 
emission profile would be fairly similar to the absorption profile. 

5. Collisional Excitation 

For the N ~ 10"^^ cm~^ clouds that we have been considering, cloud cooling is 
likely to be dominated by emission from coUisionally excited HI (see e.g., figs. I 
and 2 of Katz, Weinberg, & Hernquist 1996). Here we show that for such clouds, 
~ 50% of the energy of ionizing photons is emitted in the form Lya photons. 
For simplicity of exposition, we initially consider that the ionizing spectrum is 
restricted to < p < Apq. We include photons with p > Apq below and show that 
~ 50% of their energy also emerges in Lya. 

Consider an ionization by a photon with pq < p < Apq. The resulting electron 
has energy h[p — pq)- The electron thermalizes and (assuming ionization equilib- 
rium) is eventually recaptured when it has kinetic energy hpx ^ kT. Hence, the 
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original energy of the photon can be divided into two pieces: the energy (relative 
to the ground state) of the recombination electron, hPr^ and the thermal energy 
donated to the gas, hpg. Explicitly, 

Pr = 1^0 + VT-, Vg = V - VQ - VT- (5.1) 

The final results depend only weakly on the mean kinetic energy at recombination. 
We henceforth adopt (z/y) = O.Iz/q. 

Our basic line of argument is that about half of each of these two pieces of en- 
ergy is separately converted into Lya photons. We focus first on the recombination 
piece. Recall that each ionization leads to ~ 0.62 Lya photons. Hence the fraction 
of the recombination energy which emerges in Lya is ~ 0.62 X 0.75z/o/ ('^r) = 42%. 

We now turn to the thermal energy, hvg. We assume that 10% of this energy 
is emitted in Bremsstrahlung and that the rest goes into coUisional excitation or 
ionization. We consider four classes of collisions, with final states 1) ionization, 2) 
n = 2, 3) n > 2, and 4) n ^ 3. 

If the atom is ionized, then the fraction of the energy converted into Lya is 
the same as the recombination fraction computed above, 42%. Taking account of 
the Bremsstrahlung radiation reduces this to 0.9 X 0.42 = 38%. 

For all line excitations, 10% of the energy is lost to Bremsstrahlung, and 4% of 
the resulting Lya photons are lost to dust (just as is true of Lya photons generated 
by recombination). For n = 2, 25% go to 25* and 75% to 2P. Hence the energy 
fraction in Lya is 0.9 X 0.96 X 0.75 = 65%. For n ^ 2, excitations are mostly 
into high angular momentum states. Hence, the penultimate step of the cascade 
is nearly always 2P, which always yields a Lya photon. However, ~ 25% of the 
energy is lost in the cascade, so again the net fraction of energy in Lya is 65%. 

Finally, consider n ^ 3. For n = 3, a fraction 5/32 of the energy is lost to 
Balmer photons, and only 2/3 of the resulting n = 2 states are 2P. Hence, 9/16 of 
the excitation energy goes to Lya for a net fraction 49%. For n > 3, the fraction 
rises until it eventually reaches the n ^ 2 limit of 65%. 
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In fact, for plausible temperatures T < 30, 000 K, there are very few ionizations, 
and most excitations are to n = 2. We therefore estimate that ~ 60% of hvg is 
emitted in Lya. With this estimate we find that the net fraction of ionization 
energy emitted in Lya is 

?7£ ~ 0.50 + 0.1(^1-^^, (5.2) 

where h (^v) is the mean ionizing photon energy. 

We now take account of energetic photons, with v > Apq. These require two 
types of correction, but both are small for the HM spectrum. First, if the photons 
are sufficiently energetic, then they will pass through the cloud without ionizing 
an atom. We identify the energy at which this is a significant effect by comparing 
the optical depth to Hell of a high energy photon with the optical depth to HI 
of a photon at the Hell ionization threshold, 4z/o- The two are equal at ;/ ~ 
16(r/4)^/"^"^ z/Q ~ 40z/o where r is the ratio of Hell to HI and where the evaluation 
is for r = 50. Thus, if a cloud is optically thick for p < 4z/o, it will also be optically 
thick for p < AOpq. The precise value of this limit is not very important for the 
HM spectrum, since the absorption by the Lya forest leaves relatively few photons 
above 4z/o in any case. 

A second effect is that for each Hell ionization, a net of ~ hpo is lost during the 
process of Hell recombination. For the HM spectrum, only ~ 5% of the ionizing 
photons have p > 4z/o, while the mean frequency is {p) ~ 2.37z/o- Hence, Hell 
recombination causes a loss of ~ 2% in the energy available for Lya emission. We 
therefore ignore this effect. 

If the cloud temperature is sufficiently low, then free-free and free-bound cool- 
ing may compete with line cooling, but this should only happen if the thermal 
energy input from the ionizing background is very low, in which case the Lya 
emission from recombination would dominate over that from excitation anyway. 
Line cooling also becomes inefficient at high temperatures, but a hot cloud would 
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be coUisionally ionized and unlikely to be an optically thick system in the first 
place. Thus the conclusion that ~ 50% of the incident energy of the ionizing 
background emerges in the form of Lya photons should be fairly robust. 

We have assumed that the only energy input to the cloud is from photoion- 
ization and that other heating mechanisms, such as shocks dissipating internal 
motions, are unimportant. If other mechanisms do heat the cloud, some of this 
energy will also emerge in line cooling, enhancing the Lya emission. 

6. Observing the Forest 

6.1. Signal and Noise 

Combining the above argument for the total emission rate with the discussion 
in § 2 of the observed photon fiux, we find 



Here rj^ is the fraction of energy absorbed from the ionizing background that is 
emitted in the form of Lya photons. As argued in § 5, one expects ?7£ ~ 0.5, but 
given the range of plausible physical conditions, rj^ could be up to 20% higher or 
lower than this value. 

Suppose that a Lya cloud of angular extent AO is observed for a time T by 
a telescope of aperture (7r/4)_D^ and efficiency /. The total number of photons 
collected will then be Hs = (7r/4)C/£'^rAO$obs, where ( ~ 80% is the fraction of 
photons penetrating the Earth's atmosphere. To determine the number of noise 
photons, n„, we assume (based on the profile calculation in § 3) that the effective 
width of the double-peaked emission line is [J dv[g{v)] 

21-1 = 47ri/V , where a is the 
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1- dimensional cloud velocity dispersion and g[v) is the profile shape. We assume 
that the spectrum is taken in the B band with adequate resolution to sample the 
velocity structure of the line and with sky of i? = 22.2 and hence a sky flux of 
^sky = 1-8 X f O~^photons m~^ s~"^ A~"^ arcsec"^. The number of sky photons is 
then n„ = (7r/4)_D^/rAO^skyA/c J dv[g{v)]'^, where A is the observed wavelength. 
Hence, the signal-to-noise ratio for an optically thick cloud is, 



s . $ 



obs 



N "(1 + ^)1/2 



16(cr/c)^skyALyc 



1/2 



-3.5 / ^ \ l/2 / n \ / T \l/2 / AO \l/2/ ^ x-1/2 



l + z\ f f \ ' f D \ f T Y'' f AO \ ' / a 



3.2 J \0.25 J VlOmy V20hourJ VlOarcsecV V35km 

(6.2) 

Note that the evaluation is for a twin-peaked profile characteristic of a homogeneous 
cloud. For a Gaussian profile (which might be generated by an ensemble of small 
clouds traveling at high Mach number - see § 5) the same S/N is achieved for twice 
the dispersion, i.e. a = 70kms~"^. 

In making our numerical evaluation, we have adopted the values {p) = 2.37po 
and dv<f){v) = 5.3 X 10^ cm~^ s~"^ sr~"^ based on the HM spectrum. One would 
obtain the same values from a background with energy spectrum J{i') = 6.1 X 
10~^^ (zz/z/q)""^'^^ erg cm~^ s~"^ sr~"^ Hz~"^ . Uncertainties in the amplitude of the 
ionizing background introduce a factor ~ 2 uncertainty in the predicted S/N. 
In particular, if star-forming galaxies or dust-obscured quasars contribute to the 
background, the signal could be significantly higher. From the form of equation 
(6.2), it is clear that the prospect of detecting Lya emission declines rapidly with 
redshift, and therefore that observations should first be attempted at the minimum 
redshift where there is good system throughput. 
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6.2. Bare Detection 



As indicated by equation (6.2), even to detect Lya emission requires a long 
integration of an optically thick cloud with a high-throughput spectrograph on 
a large telescope. For an initial experiment, it is essential therefore to choose a 
system which is known to exist and to be optically thick from its absorption of 
quasar light. At ~ 10"^^cm~^, ~ galaxies observed at 2; ~ 0.5-1.5 produce 
absorption out to ~ 40/i~"^ kpc spherical radius (Steidel & Dickinson 1995), where 
h is the Hubble parameter in units of lOOkms""^ Mpc~"^. Since 1" kpc (at 

z=2.2), this corresponds to a projected diameter ~ 20". Simulations (Katz et al. 
1995) suggest that at A^ ~ 10"^^ cm~^ the clouds are smaller by a factor ~ 2. Thus 
a long 1" slit might cover ~ lOarcsec^ of a cloud. 

For an initial detection one should therefore place as wide a slit as possible over 
a quasar with a known system at z ~ 2.2 along the line of sight. The cloud will 
itself absorb most of the quasar light near Lya, but to minimize contamination, 
one should choose the quasar to be as faint as possible. If the system does extend 
over many square arcsec, it should contain emitting regions well away from the 
quasar. 

To maximize the S/N, one would like to resolve the lines. In fact, just to 
match the expected width of the lines a ~ 35kms~"^ (as assumed in making the 
S/N estimates in § 3) requires a resolution ~ 8,000. Still higher resolution would 
be better, but it is probably inconsistent with high throughput at the present time. 

A measurement of the Lya surface brightness from an optically thick cloud 
would yield a direct estimate of the integrated energy flux of the ionizing back- 
ground through equation (6.1). As shown in § 3, the emitted surface brightness is 
virtually independent of cloud geometry, so this estimate should be robust. The 
principle theoretical uncertainty is the possible contribution of heating sources 
other than photoionization; if other sources are important, then the estimate of 
the background energy density becomes an upper limit. This issue can be in- 
vestigated empirically by comparing results from a number of systems. If they 
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consistently yield the same surface brightness, then it is likely that internal heat- 
ing is unimportant and one has estimated the background energy density. If they 
yield a range of values, then one has an upper limit on the background and some 
information about the magnitude and range of internal heating sources. 

6.3. Mapping High-Redshift Structure 

At present the Lya forest is generally studied in only one dimension: along 
isolated lines of sight to quasars. These lines typically contain of order 2-3 Lyman 
limit system per unit redshift. If the clouds are indeed 5" in radius, then a l."5 X 
1000" slit should contain 500 systems per unit redshift. Assuming 2 pixels per 
resolution element, a 2048^ CCD could cover a redshift range A2; ~ 0.4 near 
z ~ 2.2. Hence, one could hope to find 200 new systems as part of the same 
observations used to detect a known absorption system in emission. This would 
provide a good statistical measure of the number density and sizes of such systems. 

In addition, once the possibility of finding such field systems were demon- 
strated, one could initiate searches in regions of the sky not associated with quasars. 
By studying the distribution of Lya clouds in two dimensions (redshift and posi- 
tion along slit), one could obtain important constraints on their correlations and 
geometry, and ultimately on their physical nature. A more ambitious project would 
be to take spectra along adjacent slit positions and so build up a 3-dimensional 
picture of the Lya forest. 

Observations like these could also "hitchhike" on programs requiring long dura- 
tion, high throughput spectroscopic exposures of the appropriate resolution. They 
could easily go hand-in-hand with searches for primeval galaxies. Indeed, it is 
quite possible that the first large population of objects detected in such searches 
will not be galaxies emitting the energy produced by their own stars but hydrogen 
gas clouds fiuorescing in the UV background. 

Acknowledgements: We thank Jordi Miralda-Escude for pointing out the 
important role of line cooling emission and for helpful discussions about possible 



21 



topologies of the emitting regions. We thank James Lowenthal and Craig Hogan 
for helpful discussions and references to earlier work. We thank Piero Madau for 
providing the Haardt & Madau (1995) spectrum in a convenient numerical form, 
prior to publication. Work by A.G. was supported in part by grant AST 94-20746 
from the NSF and by grant NAG5-3111 from NASA. Work by D.W. was supported 
in part by NASA grants NAG5-2882 and NAG5-3111. 



22 



REFERENCES 

Fall, S. M. & Pei, Y. C. 1995, in QSO Absorption Lines, Proc. ESO Workshop, 
ed. G. Meylan (Heidelberg: Springer), p. 23 

Haardt, F. & Madau, P. 1995, ApJ, in press (HM) 

Hogan, C. J. & Weymann, R. J. 1987, MNRAS, 225, PI 

Holstein, T. 1947, Pliys. Rev. 72, 1212 

Katz, N., Weinberg, D. H., Hernquist, L., & Miralda-Fscude, J. 1995, ApJ Letters, 
in press 

Katz, N., Weinberg, D. H., & Hernquist, L. 1996, ApJ Supp, in press 

Lowentlial, J., Hogan, C. J., Green, R. F., Woodgate, B., Caulet, A., Brown, L., 
& Bechtold, J. 1995, ApJ, 451, 484 

Lowenthal, J., Hogan, C. J., Leach, R. W., Schmidt, G. D., & Foltz, C. B. 1990, 
ApJ, 357, 3 

Mo, H. 1995, MNRAS, in press 

Neufeld, D. A. 1991, ApJ, 370, L85 

Osterbrock, D. F. 1989, Astrophysics of Gaseous Nebulae and Active Galactic 
Nuclei (Mill Valley, University Science Books) 

Steidel, C. C. & Dickinson, M. 1995, in Wide Field Spectroscopy and the Distant 
Universe, eds. S. J. Maddox & A. Ar agon- Salamanca (Singapore: World 
Scientific) 

Thompson, D., Djorgovski, S., & Trauger, J. 1995, AJ, 110, 963 
Thompson, D. & Djorgovski, S. 1995, AJ, 110, 982 



23 



